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Abstract 

Consider the optimal stopping problem of a one-dimensional diffusion with posit- 
ive discount. Based on Dynkin's characterization of the value as the minimal excessive 
majorant of the reward and considering its Riesz representation, we give an explicit 
equation to find the optimal stopping threshold for problems with one-sided stop- 
ping regions, and an explicit formula for the value function of the problem. This 
representation also gives light on the validity of the smooth fit principle. The results 
are illustrated by solving some classical problems, and also through the solution of: 
optimal stopping of the skew Brownian motion, and optimal stopping of the sticky 
Brownian motion, including cases in which the smooth fit principle fails. 



1 Introduction and problem formulation 

Consider a non-terminating and regular one-dimcnsional (or linear) diffusion X = {Xt : t > 
0}, in the sense of Ito and McKean [9 (see also Borodin and Salminen [2 ). The state space 
of X is denoted by I, an interval of the real line M with left endpoint £ — mil and right 
endpoint r = sup I, where — oo < £ < r < oo. We exclude the possibility of absorbing and 
killing boundaries; if some of the boundaries belong to I we assume it to be both entrance 
and exit (i.e. non-singular). Denote by P^; the probability measure associated with X 
when starting from x, and by the corresponding mathematical expectation. Denote by 
T the set of all stopping times with respect to {J-t ■ t > 0}, the usual augmentation of the 
natural filtration generated by X (see 1.14 in [?:). 

Given a non-negative lower semicontinuous reward function g: I ^ R and a discount 
factor a > 0, consider the optimal stopping problem consisting in finding a function 
and a stopping time r* G T, such that 

Va^ix) = (e-""*.g(X,.)) = sup (e-"^*g(X,.)) • (1) 

The value function Va and the optimal stopping time t* are the solution of the problem. 

The first problems in optimal stopping appeared in the framework of statistics, more 
precisely, in the context of sequential analysis (see the book by Wald [29 ) . For continuous 
time processes a relevant reference is the book of Shiryaev [27] that also has applications 
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to statistics. A new impulse to these problems is related to mathematical finance, where 
arbitrage considerations give that in order to price an American option one has to solve 
an optimal stopping problem. The first results in this direction were provided by Mc Kean 
[H] in 1965 and Merton [TS] in 1973, who respectively solved the perpetual put and call 
option pricing problem, by solving the corresponding optimal stopping problems in the 
context of the Black and Scholes model [1]. For the state of the art in the subject see the 
book by Peskir and Shiryaev [20^ and the references therein. A new approach for solving 
onc-dimensional optimal stopping problems for very general reward functions is provided 
in the recent paper |19j . 

When considering optimal stopping problems we typically find two classes of results. 
The first one consists in the explicit solution to a concrete optimal stopping problem ([I]). 
Usually in this case one has to -somehow- guess the solution and prove that this guess 
in fact solves the optimization problem; we call this approach "verification" . For example 
we can consider the papers [13], [H], [IH], [li]. The second class consists of general 
results, for wide classes of processes and reward functions. We call this the "theoretical" 
approach. It typically include results about properties of the solution. In this class we 
mention for example [3], [8], [6]. But these two classes not always meet, as frequently in 
concrete problems the assumptions of the theoretical studies are not fulfilled, and, what 
is more important, many of these theoretical studies do not provide concrete ways to find 
solutions. In what concerns the first approach, a usual procedure is to apply the principle 
of smooth fit, that generally leads to the solution of two equations: the continuous fit 
equation and the smooth fit equation. Once these equations are solved, a verification 
procedure is needed in order to prove that the candidate is the effective solution of the 
problem (see chapter IV in .20j). This approach, when an explicit solution can be found, is 
very effective. In what concerns the second approach, maybe the most important result is 
Dynkin's characterization of the solution of the value function Va as the least a-excessive 
(or a-superharmonic) majorant of the payoff function g [3 . Other ways of classifying 
approaches in order to study optimal stopping problems include the martingale-Markovian 
dichotomy as exposed in [5D]. 

In the present work we provide an explicit solution of a right-sided optimal stopping 
problem for a one- dimensional diffusion process, i.e., when the optimal stopping time has 
the form 

T* ^mf{t>0: Xt>x*} (2) 

for some optimal threshold x* G I, under mild regularity conditions. Right-sided problems 
are also known as call-type optimal stopping problems. Analogous results are valid for left- 
sided problems. 

An important byproduct of our results has to do with the smooth fit principle. Our 
results are independent of this principle, but they give sufficient conditions in order to 
guarantee it. 

In section [2] some necessary definitions and preliminary results are given. Our main 
results are presented in section |3] In section |4] we discuss the consequences of our results 
related to the smooth fit principle. Finally, in section [5] we present some examples, includ- 
ing the optimal stopping of the skew Brownian motion and of the sticky Brownian motion 
(suggested by Paavo Salminen), where particular attention to the smooth fit principle is 
given. 
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2 Definitions and preliminary results 



Denote by L the infinitesimal generator of the diffusion X, and by 2?^ its domain. For 
any stopping time t and for any / G 2?^ the following discounted version of the Dynkin's 
formula holds: 

f{x) = (^^" e-"*(« - L)f{Xt)dt^ + E, (e-"V(^r)) • (3) 



The resolvent of the process X is the operator defined by 

R^u{x)= / e-"*E, 
Jo 



applied to a function u G Cf,(I) = {u: I — >■ K, m is continuous and bounded}. The range of 
the operator Ra is independent of a > and coincides with the domain of the infinitesimal 
generator V^. Moreover, for any / e Ra{a — L)f = /, and for any u <E Cb(I), 
(a — L)RaU = u. In other terms, Ra and (a — L) are inverse operators. Denoting by s 
and m the scale function and the speed measure of the diffusion X respectively, we have 
that, for any f £ Vl, the lateral derivatives with respect to the scale function exist for 
every x £ {(,r). Furthermore, they satisfy 

^ix)-^ix)=m{{x})Lfix), (4) 
and the following identity holds for z > y: 



^^^-{z)-K^{y)^ I Lf{x)m{dx). (5) 



ds ds 



This last formula allows us to compute the infinitesimal generator of / at a:; S (i, r) by the 
Feller's differential operator [7] 

The infinitesimal generator at £ and r (if they belong to I) can be computed as Lf{£) = 
lim^_^i+ Lf{x) and Lf{r) = lim^. Lf{x) respectively. 

Given a function u: X — > M, and x G {i,r) we give to Lu{x) the meaning given in 
^ if it makes sense. We also define, if ^ G X, Lu{i) = lim2,_j.£+ Lu{x) and if r e I, 
Lu{r) = lim^_j.r~ Lu{x), if the limit exists. 

There exist two continuous functions Lpa'- 1 ^ IR+ decreasing, and '0q '■ ^ hi- 
creasing, solutions of au — Lu, such that any other continuous function u is a solution 
of the differential equation if and only if tt = aipa + hif^on with a,b in M. Denoting by 
Tz = infjt : Xt = z} the hitting time of level z e I, we have 




The functions and tpa, though not necessarily in V^, also satisfy Q for all x G 
which allow us to conclude that in case m{{x}) = 0, the derivative at x of both fmictions 
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with respect to the scale exists. From ([5| apphed to tpa, and taking into account af/'a 
Lip a we obtain for z < y 

"(^) ^~^^y)= / atpa{x)m{dx); 



ds ds 



the right hand side is strictly positive since aipa is positive and m charge every open 
set. We conclude that — gj^ is strictly increasing. In an analogous way it can be proven 
that is increasing as well. The previous consideration, together with the fact that 

^ (x) > and (x) < allow us to conclude for x G {i,r): 

-oo < — — (x) < — — (x) < < — (x) < — (x) < oo. 

OS OS OS OS 

The Green function of the process X with discount factor a is defined by 

/•oo 

Ga{x,y) = e''°'*p{t;x,y)dt, 



where p{t; a;, y) is the transition density of the diffusion with respect to the speed measure 
m{dx) (this density always exists, see [9] or [2]). The Green function may be expressed in 
terms of ipa and ipa as follows: 

G (x y) = X <y, 

\w~'^1pa{y)^a{x), X>y, 

where Wa is the Wronskian, given by 

Wa = ^^^{x)(pa{x) - ^a{x)^^{x). 

Observe that the Wronskian is positive and independent of a; Given u: I R, 

under the condition J^- Gaix,y)\u{y)\'m{dy) < oo, an application of Fubini's Theorem 
gives 

Rau{x) ^ J Gaix,y)u{y)m{dy). (9) 

A non-negative Borel function u: I — > M is called a-excessive for the process X if 
e~"* {u{Xt)) < u{x) for all x e 1 and t > 0, and limt_j.o E^; {u{Xt)) — u[x) for all 
a; € I. A 0-excessive function is said to be excessive. 

Consider the process killed at an independent exponential time of parameter a, i.e. 




with Ca a random variable with exponential distribution of parameter a independent of 
X, and A the cemetery state, at which any function is defined to be zero. It is easy to 
see that the Green function G^ of the process Y — {Yt : t > 0} coincides with Gq,; a Borel 
function u: T i— ^ M is excessive for Y if and only if it is a-excessive for X. In fact, the 
non-discounted optimal stopping problem for the process Y has the very same solution 
(value function and optimal stopping time) as the a-discounted optimal stopping problem 
for X. 

For general reference on diffusions and Markov processes see [2 151 [211 El HH] • 
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3 Main results 



Our departing point, inscribed in the Markovian approach, is Dynkin's characterization 
of the optimal stopping problem solution. Dynkin's characterization M states that, if the 
reward function is lower semi-continuous, V is the value function of the non-discounted 
optimal stopping problem with reward g if and only if V is the least excessive function 
such that V{x) > g{x) for all x £ I. Applying this result for the killed process Y, and 
taking into account the relation between X and Y, we obtain that V^, the value function 
of the problem with discount a, is characterized as the least a-excessive majorant of g. 

The second step uses Riesz's decomposition of an a-excessive function. We recall this 
decomposition in our context (see [HI [12 HI)- A function u: Z ^ M is a-excessive if and 
only if there exist a non-negative Radon measure and an a-harmonic function such that 

u{x) = / 6*0,(2;, y)/i(c??;) -1- (a-harmonic function). (10) 
J{e.r) 

Furthermore, the previous representation is unique. The measure /i is called the repres- 
enting measure of u. 

The third step is based on the fact that the resolvent and the infinitesimal generator 
of a Markov process are inverse operators. Suppose that we can write 

V^{x) = j^G^{x,y){a~ L)V^{y)m{dv), (11) 

where L is the infinitesimal generator and m{dy) is the speed measure of the diffusion. 
Assuming that the stopping region has the form T f] {x > x*}, and taking into account 
that Va is a-harmonic in the continuation region and Va = g the stopping region we 
obtain as a suitable candidate to be the representing measure 

{0, ify<2:*, 
M,.(d2/), iiy^x\ (12) 

{a ~ L)g{y)m{dy), if y > x* , 



This approach was initiated by Salminen in [52] (see also [Hj). According to Salminen's 
approach, once the excessive function is represented as an integral with respect to the 
Martin kernel Ma{x,y), 

Vc{x) ^ J^M^{x,y)K{dy) (13) 
one has to find the representing measure k. Martin and Green kernels are related by 

Ga{xo,y) ■ 

Ga(xo,V) 



Ma{x,y) = (?°/^'^\ ; where xq is a reference point. Therefore, the Riesz's representation 



of Va is related with the one in ( 13 ) by considering iy{dy) = g'^\x^\) Ma {x, i)K{{£}) ■ 



Maix,r)K{{r}) as the a-harmonic function. 

It is useful to observe that when the optimal stopping problem ([l]) is right-sided with 
optimal threshold x* it has a value function Va of the form 



Vaix) 



Ex (e""""-*) 5(2;*), x<x* 
g{x), x>x* 
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Furthermore, Va{x) > g(x) for all x € T and, in virtue of equation ([7|), we have 



9{x), 



X > X 



(14) 



The state space of the process can include or not the left endpoint £ and the right 
endpoint r. In order to simplify, with a slight abuse of notation, we write [i,x], [£,x), 
[x,r], {x,r] to denote respectively I Ci {y < x}, X Ci {y < x}, In {y > x}, X Ci {y > x}. 

We say that the function g: X R satisfies the right regularity condition (RRC) if 
there exist a point Xi d X and a function 5 : I — > M (not need to be non-negative) such 
that g{x) — g{x) for x > xi and 



g(x) = / Ga{x,y){a- L)g{y)m{dy) {xeX). 



(15) 



Proposition 3.5 gives conditions in order to verify the inversion formula (15 1. Informally 
speaking, the RRC is fulfilled by functions g that satisfy all the local conditions -regularity 
conditions- to belong to I?l for x > xi, and does not increase as quick as tpa does when 
approaching r (in the case r ^ X). Observe that if g satisfies the RRC for certain xi it 
also satisfies it for any greater value; and of course, if g itself satisfy ( 15 1 then it satisfies 
the RRC for all xi in X. To take full advantage of the following result it is desirable to 
find the least xi such that the RRC holds. 

The main result follows. 
Theorem 3.1. Consider a diffusion X and a reward function g that satisfies the RRC 
for some xi. The optimal stopping problem is right-sided with optimal threshold x* > xi 
if and only if: 



g{x*) > w^'^ipaix*) 



Va{y){a - L)g{y)m{dy), 



and 



' (x* .r] 

{a — L)g{x) > 0, x £X: x > x* 
](.x^ ) 

-ifja{x) > g{x), X GX: X < X* 



(16) 
(17) 
(18) 



Tpaix*) 

Furthermore, in the previous situation: 

• the Riesz 's representation of the value function Va has representing measure as given 



in ( 12 I with 



g{x*) - Va(a;*)/(^._^] <Pa{y){a ~ L)g{y)m{dy) 



Va(a;*)<y3Q(x*) 

while the a-harmonic part vanishes; 



• if X* > Xi and the inequality (16) is strict, then x* is the minimal satisfying this 
strict inequality and (17), in particular 



g{x*) < ^-ipaix*) / Vaiy)ia - L)g{y)m,{dy); 



X* ,r] 



(19) 



which implies that k < (a — L)g{x*)m{{x*}); 
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Remark 3.2. From this theorem we obtain an algorithm to solve right-sided optimal 
stopping problems which works in most cases: (i) Find the largest root x* of the equation 



g{x*) = w^'^i^aix*) I Lpa{y){a - L)g{y)m{dy); 



(20) 



(ii) Verify (a — L)g{y) > for x > cc*; (iii) Verify g{x) < Jj^^^l-j 'ipaix)- If these steps are 
fulfilled, the problem is right-sided with optimal threshold x* . Observe that if m{{x*}) 
0, then inequalities ( 16 1 and ( 19 ) are equalities; 



Proof. We start by observing that if the problem is right-sided with threshold x* then (17) 



holds. In general {a — L)g is non-negative in the stopping region (this can be seen with 
the help of the Dynkin's operator, see ex. 3.17 p. 310 in [21], see also equation (10.1.35) 
in [13 )■ Under the made assumption the value function Va is given by (14 1, being ( [Ts] ) a 
consequence of the fact that the value function dominates the reward. To finish the proof 



of the "only- if" part it remains to prove (16 1. Consider Wa : 1 1— >■ M defined by 



W^ix) 



{x* ,r] 



Ga{x, y){a - L)g{y)m{dy); 



observe that is a-excessive in virtue of and Riesz's representation. Let : X i-> M 
be defined by 

Vaix) := Wo^ix) + kGa{x,x*), 

where k is such that Va{x*) = g{x*), i.e. k = {g{x*) - Wa{x*)) /Ga{x* ,x*). Observe 
that, by Wa{x*) is the right-hand side of ([l6|. In fact, ^ holds if and only if fc > 0. 
By the definition of Va and the representation (|8|) of Ga we get for x < x* 

Let us compute Va{x) — g{x) for x > x* . In this region we have g = g, where g is 



the extension given by the RRC. For g can use the inversion formula (15). Denoting by 

Vh{dy) = (a — L)g{y)m{dy) we have 

Va{x) - g{x) ^ kGa{x,x*) - / Ga{x,y)vh{dy) 

J[l,x*] 

iPa{x) 



(fiaix*) 



{Va{x*)~gix*)) = 0, 



because Va{x*) = g{x*). So far, we have proved that Va{x) = Va{x) for all x <El. We 
are ready to prove that fc > 0, based on the uniqueness of the Riesz's decomposition: the 
a-excessive function Wa has the Riesz representation given by its definition, and, if fc < 
then 

Wa{x) = -kGa{x,X*)+Va{x) 

would give another Riesz's representation (the representing measure being —k5^^*^{dx) -\- 
fj,{dx), where fj, is the representing measure of Vq). An easy way of verifying that both 
measures are not the same is to observe that the former docs not charge {x*}, while the 
latter do. 



7 



To prove the "if" statement observe that, assuming (16) (17) and (18), function Va, 
already defined, is a-excessive (by Riezs's representation, bearing in mind that fc > 0) 
and dominates g. By Dynkin's characterization, the value function Va is the minimal 
a-excessive function that dominates g. Therefore Va > Va- Since Va satisfies (14) (we 
have proved this in the fist part of this proof), it follows that Va is the expected reward 
associated with the hitting time of the set [x* , r], then 

Vaix) < Vaix) = inf E, (e-"^5(X,)) , 

r 

concluding that the problem is right-sided with threshold x* . 

The consideration about the Riesz's representation of Va stated in the "furthermore" 
part are a direct consequence of the made proof. 

To prove the minimality of x* , suppose that there exists x** such that xi < x** < x* 
satisfying the strict inequality in ([l6| and ([l7|. Let us check Va{x**) — g{x**) < 0, in 
contradiction with the fact that Va is a majorant of g. Considering the extension g given 
by the RRC and denoting Vg{dy) — [a — L)g(jj)m{dy) we have 

Va(x**) ~ g{x**) = - / Ga(x*\y)Vg{dy) + kGa{x*\x*) = Si+S2 + Sj, + Si, 



where 



si = - / Ga{x** ,y)v~g{dy), 

J[l,x"] 

32 ^ -Wa^Tpaix**) Lp a{y)Vci{dy) , 

J{x-",x*] 



Ipaix*) iPa{x**) J[i,^*,] 



*4 = t^*-. / Ga{x*,y)i^g{dy). 

YayX ) J(x",x''] 

To check that S3 -I- S4 = kGa(x** ,x*), use 

, 1 [ n ( * ^ M \ A Ga{x*\x*) i;a{x**) 

^ = TTTZTZIK f Ga{x ,y)vg{dy), and 



Ga{.X* ,X*) Jyn^^.} Ga{x*,X*) tpa{x*) 

Finally, observe that 



i^ajx**) (fiajx*) 
i'aix*) (Pa{x**) J[e^x"] 



^i+-^3 = ( r7,.^:7:J^ -l) / G„(x-,y)^,(dy)<o, 



because the first factor is negative and the second one positive, by the assumption about 
a;**, and 



W g^lpajx**) 

4'a{x*) J (x" .x*] 



*4 + S2 = °,, / {tp{y)ipa{x*) - ^a{x*)ipa{y)) V~g{dy) < 0, 



because the measure Vg (dy) is positive and the integrand non-positive (it is increasing and 
vanishes in y — x*). We have obtained that Va{x**) — g{x**) = si + S2 + S3 + S4 < 0, 
concluding the proof. □ 



8 



The previous theorem gives necessary and sufRcient conditions for the problem ([I]) to 
be right-sided under the RRC. The following result gives simpler sufficient conditions for 
the problem to be right-sided. 

Theorem 3.3. Consider a diffusion X and a reward function g such that (15 I is fulfilled 
for g = g. Suppose that there exists a root c £ f ) of the equation {a — L)g{x) = 0, 
such that (a — L)g{x) < if x < c and {a — L)g{x) > if x > c, and that Jj-'ipa{y){a — 
L)g[y)m{dy) G (0, oo]. Then the optimal stopping problem ([l]) is right-sided, with optimal 
threshold 

X* ^uAn{x:h{x)>Q), (21) 



whe 



h{x) 



ijja{y){a- L)g{y)m{dy) [x el). 



Proof. The idea is to apply Theorem 3.1 with x* defined in (21 ). By the assumptions on 
(a — L)g and the fact that m{dy) is strictly positive in any open set, we obtain that the 
function h{x) is decreasing in [£, c) and increasing in (c, r). Moreover h{x) < Q \i t < x < c. 
Since 6 is right continuous and increasing in (c, r), the set {x: b{x) > 0} = [x* ,r) with 
x* > a. Observe that, by (15) and (|8| we get 



g(x) = W^'^ipaix)b{x) 



{x,r] 



Ga{x,y){a - L)g{y)m{dy), 



and b{x*) > is equivalent to (16 1. Since a;* > c we have {a — L)g{y) > for x > a;*. It 
only remains to verify ( 18 1. By definition of x* there exists a signed measure at{dy) whose 
support is contained in x*], and ai(dy) — {a — L)g{y)m(dy) for y < x* and such that 



ipaiy)(Te{dy) = 0. 



Furthermore ar{dy) — {a — L)g{y)m{dy) — a p{dy) is a positive measure supported in [a;*, r] 
Using the inversion formula for g and (18]), we have for x < x* 



5W - , / ^^ 9{x ) 

VaKX*) 



G^ix,y)a,idy)<^^^ 



i^a{y)(^i{dy) = 0, 



where the inequality follows from the following facts: if y < c then ai{dy) < and 



^a(y)%^7T^ < Ga{x,y), 



while if 2/ > c then (T£{dy) > 



^a(y)^#^>G„(x,y). 



i^a{c) 

We can now apply Theorem |3 . 1 1 completing the proof. 



□ 



3.1 On the right regularity condition (RRC) 



In order to apply the previous results it is necessary to verify the inversion formula (15 1. 
As we have seen in the preliminaries, if / G we have Ra{a — L)f =^ /, and if equation 
([9| holds for {a — L)f, we have (15). This is the content of the following result. 
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Lemma 3.4. Assume f G Vl, and 

Gai^, y)\(,ct ~ L)f{y)\m{dy) < oo for all x £l. 



Then equation (15) holds. 

The conditions of the previous lemma are very restrictive in order to solve concrete 
problems, as reward functions typically satisfy lim^^-j.^ 5(2:) = 00. The following result 
extends the previous one to unbounded functions. 

Proposition 3.5. Consider the case r ^ I. Suppose it: 1 M is such that Lu{x) in ^ 
can he defined for all x € I. Assume 

Gaix,y)\{a~ L)u{y)\m{dy) < 00, (22) 

Jx 

and 

hm ^ = 0. (23) 

Suppose also that for each y G I there exist a function Uy G I?l such that Uy{x) — u{x) 
for X < y. Then u satisfies (15). 

Proof. By ([9]) we have Jj-Ga{x,y){a~ L)u{y)m{dy) — Ra{a — L)u{x). Consider a strictly 
increasing sequence ?■„ — > r (n — > co) and denote by u„ the function Ur,,^^ G I?l of 
the hypothesis. By the continuity of the sample paths, by our assumptions on the right 
boundary r, we have Tr„ — >■ 00 (ri — >■ 00). Applying formula ([3| to u„ and the stopping 
time Tr,^ we obtain, for a; < r„, 

Unix) = E, ^ e'"*(a - L)un{Xt)dt^ + E, (e'""-" ) M„(r„), 

using Unix) — uix) and (a — L)u{x) — ia — L)w„(a;) for x < r„+i we have 

uix) = E, ^ e-"*(a - L)uiXt)dt^ + E, (e'""- ) u(r„). 

Taking limits when n ^ cx), by ([t]) and (23) we have 

E. (e-"^-) i.(r„) = ^I^Mrn) ^ 0. 

To compute the limit of the first term above we use dominated convergence theorem and 
( [22| ). The result is 

u{x)= / E, (e-"*(a-L)u(XO)rft = / G„(a;, y)(a - i)u(y)m(dy) 



concluding the proof. □ 



4 On the principle of smooth fit 

The principle of smooth fit (SF) holds when condition V'ix*) = g'ix*) is satisfied, being 
a helpful tool to find candidate solutions to optimal stopping problems. In 22] Salminen 
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proposes an alternative version of this principle, considering derivatives with respect to 
the scale function. We say that there is scale smooth fit (SSF) when the value function 
has derivative at x* with respect to the scale function. Note that if g also has derivative 
with respect to the scale function they coincide, since g = Va in [x*,r]. In [18 Peskir 
presents two interesting examples: one of them consists on the optimal stopping problem 
of a regular diffusion with a differentiable payoff function in which the principle of SF does 
not hold, but the alternative principle of SSF does; while in the other the principle of SF 
holds but the principle of SSF fails. Later, Samee |i24j analysed the validity of the principle 
of smooth fit for killed diffusions and introduced other alternative principle considering 
derivatives of and ^ with respect to the scale function s. 

We now analyse the relation between the Riez's representation of Vq, stated in the 
previous section, and the principle of smooth fit. We start by proving that k = j/({a;*}) = 
in (12 1 implies that the reward function has derivatives with respect to the function ipa- 
Then we follow by stating some corollary results. 

Theorem 4.1. Given a diffusion X and a reward function g, if the value function asso- 
ciated with the problem (jlj satisfies 



Va{x) = Gaix,y){a~ L)g{y)m{dy), 

for X* e (l,r), then Va is differentiable at x* with respect to ipa- 
Proof. For x < x* 

Va{x) =W~^1pa{x) I ipa{y)v{dy), 
(a:* ,r] 

and the left derivative of Va with respect to in x* is 

^^—{x*)=w^^[ ipa{y)v{dy). 



dtpa 

For X > X* 



(x* ,r] 



Vaix) = (pa{x)Wa^ j ti^a{y)v{dy)+^a{x)Wa^ j Lpa{y)v{dy). 

Computing the difference between Va{x) and Va{x*) we obtain 

Va{x) -Va{x*) ^ W^^ {ipaix) - i^aix*)) / (pa{y)l^{dy) 



Then 



[a;,r] 

,,-1 



/ i'^aix)i^aiy) ~ '4'»{x*)(Pa{y))l'{dy). 

J (x*- ,x) 



-1 



> + 



lim / ipa{y)v{dy) 



w„ ^ lim 



x^x* 



/(^.^^) ¥'a(a;)V'a(y) - tpa{x*)^a{y)'^{dy) 

tpa{x) - -ipaix*) 
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If the last limit vanishes, we obtain that the right derivative exists, and 



dV+ dV^ 



{x.r] 



This means that we have to prove 

f (fa {x)ll'a (y) - tpa {x*)(pa (v) 



lim 

a;— J-a;* + 



v{dy) - 0. 



(24) 



Denoting by f{y) the numerator of the integrand in (24), observe that 

.f{y) = Va{x){ll)a{y) - Tpaix*)) + ij^ a{x* ) {if a{x) - ipa{y))- 

For the first term, we have (observe that x* < y < x) 

< Lpa{x){i'a{'y) - ll^a{x*)) < ipa{x) {lpa{x) - -ipaix*)), 

while for the second 

> tpa{x*){(pa{x) - (fiaiy)) > 1pa{x*){(faix) - ipa(x*)). 

We conclude that 

ll^a{x*){(pa{x) - ipaix*)) < f (y) < ipa{x){tlJa{x) - '0q(x*)). 

Dividing by ipa{x) — ipaix*) we see that the integrand has a lower bound b{x) given by 



Ipaix) ~ 1pa{x*) 



S{x) - s{x*) Ipaix) - Ipaix*)' 



while faix) is an upper bound. We obtain the integral in (24) satisfies 

hix)vix*x)<[ Mx)i^o.iy) - 4'o.ix*)v aiy) . 

~ J(x\x) Tpaix) -Ipaix*) 

Taking limits when x — > x*^ we obtain faix) — ^ faix*), uix^x*) — > 0, and 



-vidy) < ipaix)iyix* ,x). 



lim bix)=Mx*)^ix*) /^ix*) , 

OS / OS 



concluding that (24) holds. 



□ 



As we have seen in section |3j if the speed measure does not charge x* neither does the 
representing measure. This means that if representation (12) holds, then m({a;*}) = 
is enough to guarantee the differentiability of Va with respect to ipa- We also have the 
following result. 

Corolary 4.2. Assume the conditions of Theorem \4-l\ If the speed measure does not 
charge x* there is scale smooth fit. 

Proof. By the previous theorem we know that Va is differentiable with respect to ipa- 
Condition mi{x*}) — implies that ipa has derivative with respect to the scale function. 
We conclude 



ds 



dtpa 



□ 
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The previous result, under the additional assumption that ipa and fa are differentiable 
with respect to s, could be derived from Corollary 3.7 in [22j. This result states that the 
representing measure of Va does not charg x* if and only if Va is differentiable with respect 
to s. Also Theorem 14.11 can be derived from the mentioned result under the additional 
assumption by using the chain rule. 

As a consequence of the previous results, we obtain, by using the chain rule, conditions 
under which the principle of SF holds. 

Corolary 4.3. Assume that g is differentiable at x* . Under the conditions of Theorem 



4-1: if ipa is differentiable at x* and ^/^{x*) ^ Q (or under the conditions of Corollary 4--2 



if s is differentiable at x* and s'{x*) 0) then the principle of SF holds. 

The previous result is closely related with Theorem 2.3 in [18 , which states that, in 
the non-discounted problem, there is smooth fit if the reward and the scale function are 
differentiable at x* . Theorem 2.3 in [M] ensures the validity of the smooth fit principle for 
the discounted problem under the assumption that g, ipa and ipa are differentiable at x* . 
It should be noticed that these results are valid in general, not only in one-sided problems. 



5 Examples 

In this section we show how to solve some optimal stopping problems using the previous 
results. We present two classical examples (American and Russian options), and also 
include some new examples in which the smooth fit principle is not useful to find the 
solution. 



5.1 American call options 



Consider a geometric Brownian motion given by Xt — x exp^aWt -I- (/i — tT^/2)t), where 
{Wt} is a standard Brownian motion, /i e M and cr^ > 0. The state space is / = (0,oo). 
We refer to [Ij, p. 132 for the basic characteristics of this process. The infinitesimal 



generator is Lf = icr^a;^/" -I- fJ,xf' , with domain 



VL = {f: f, Lf e Cb{I)}. 

Consider the payoff function g{x) = {x — K)~^ {x gR), where if is a positive constant, and 
a positive discount factor a satisfying a > fi. The reward function g satisfies the RRC for 
Xi = K: it is enough to consider g e C^, bounded in (0, K) and such that g{x) — x~ K for 



X > K. This function g satisfies the inversion formula ( 15 1 as a consequence of Proposition 



3.5 Observe that equation (22) holds. Equation (^23| is in this case 



lim 



Ipaiz) 



lim '^'^ , 



the last limit vanishes if 1 — 71 < 0, which is equivalent to fi < a. To find x* we solve 



equation (20). After computations, we find 



K 



71 



7i 



with 



71 




2a 
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It is not difficult to verify (171 and (18) in order to apply Theorem 3.1 We conclude 



that the problem is right-sided with optimal threshold x* . Observe that the hypotheses 
of Theorem |4 . 1 1 and corollaries |4.2| and |4.3| are fulfilled. In consequence all the variants of 
smooth fit principle hold. This problem was solved by Merton in j l5j . 



5.2 Russian Options 

The Russian Option was introduced by Shepp and Shiryaev in 1993 in |25j . where the 
option pricing problem is solved by reduction to an optimal stopping problem of a two- 
dimensional Markov process. Later, in |26) . the authors give an alternative approach to 
the same problem solving a one-dimensional optimal stopping problem. In 2000, Salminen 
[23j . making use of a generalization of Levy's theorem for Brownian motion with drift, 
shortened the derivation of the valuation formula in [26j and solved the optimal stopping 
problem related. 

Consider a > 0, r > and a > 0. Let {Xt} be a Brownian motion on / = [0,oo), 

with drift —S < 0, where S — "^^"^ and reflected at 0. In (35] it is shown that the 
optimal stopping problem to be solved has underlying process {Xt\ and reward function 
g{x) = e"^^ . For the basic characteristics of the process we refer to [2], p. 129. The 
infinitesimal generator is Lf{x) = f"{x)/2 — df'{x) for a; > and Lf{0) = limj._j.o+ Lf{x), 
with domain 

2?i = {/ : /, Lf e Chil), lim f'{x) = 0}. 

The payoff function g{x) satisfies the RRC for every xi > 0: for xi > it is easy to find a 
function g with continuous second derivatives such that g — g in [xi, oo) and such that the 



right derivative at is 0. By the application of Proposition 3.5 we obtain that g satisfies 



the inversion formula (15), then the RRC holds. We obtain 



(a - L)g{x) = {a- + Sa)e''=' = {a + r)e''^> 0. 



In order to apply Theorem 3.1 we solve equation (20) which in this case is 
1 



7 — (5 \ 27 27 
with 7 = \/2a + (5^, obtaining (observe that —7 — S + a < 0) 



X* = — In 
27 



7 + (5\ f J — 6 + a 



^ — S J \7 + (5 — (7 



It is easy to to verify conditions (17) and (18), to obtain, by application of Theorem 3.1 



that the problem is right-sided with threshold proved in |2Bj. 



5.3 Skew Brownian motion 

We consider a Brownian motion skew at zero. This process behaves like a standard 
Brownian motion outside the origin, but has an asymmetric behaviour when hitting a; = 0, 
modeling a permeable barrier. The behaviour at a; = is regulated by a parameter 
/? S (0, 1), known as the skewness parameter. The space state of this process is I = K. 
For details on this process and its basic characteristics we refer to see [2], p. 126 or [13]. 
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The infinitesimal generator is Lf{x) = f"{x)/2 ii x and Lf{0) = Mmx^o L f (x) , with 
domain 

- {/ : /, Lf e Cb(Z), /?/'(0+) = (1 - /3)/'(0-)}. 

Consider the payoff function 17(2:) — a;+. Function g satisfies the RRC for xi = 0: to see 
this it is necessary to construct g such that g — g in [0, 00) g with second derivative bounded 
in (— cx),0) and such that .g'(0~) = /3/(l ~ /3) (so that g satisfies the local conditions to 
belong to Applying Proposition |3.5| it can be concluded that g satisfies (15), so the 

RRC holds. We have (a — L)g{x) — ax {x > 0). Equation (20) is in this case 



^ ^^■sinh(V2^a;*)+e^-* 



2a V /3 



/2a t 



at 2l3dt 



or equivalently 



^ ((2/3 - l)e-^ ^* (%/2^ X* + + V2^ x* + l) 



2V2a 



(25) 



In general, this equation can not be solved analytically. If we consider the particular case 
Rrhich the process is the ordinary Brownian motion, we obtain the known result 
(see [IH])- Consider a particular case, in which a = 1 and f3 = 0.9. Solving 



1 

/2a 



numerically equation (25) we obtain x* ~ 0.82575. Checking (17) and (18) we conclude 



that the problem is right-sided with optimal threshold x* 



5.3.1 An example without smooth fit 



Consider again the Skew Brownian motion with parameter /3 — 1/3, a payoff function 
g{x) = {x + 1)+ and a discount a — 1/8. We have {a — L)g{x) = a{x + 1) {x > 0). 
Observe that x* = is a solution of (16) (with equality). It is easy to see that the 



assumptions of Theorem 3.1 are fulfilled. We conclude that the problem is right-sided 
with threshold x* = 0. Moreover, we know 



Vo^ix) 



\x + l, x>0, 

\lpa{x), X<0. 



Unlike the previous examples, the value function Va is not differentiable at x* . As we see 
in Figure [lithe graph of Va shows an angle at x = 0. By application of Theorem 4T and 
Corollary ]4. 2 1 we conclude that Va is differentiable with respect to ipa and SSF hold. 

An example considering a regular diffusion with non-differentiable scale function, in 
which the the SF fails to hold, was provided for the fist time by Peskir [TS] . 



5.4 Sticky Brownian Motion 

Consider a Brownian motion, sticky at 0. We refer to [5] p. 123 for the basic characteristics 
of this process. It behaves like a Brownian motion out of 0, but expends a positive time at 
a; = 0; this time depends on a positive parameter that we assume to be 1 in this example. 
The space state of this process is I = M. The scale function is s{x) = x and the speed 
measure is m{dx) ~ 2dx + 25^^^{dx). The infinitesimal generator is Lf{x) = /"(a;)/2 when 
X 0, and Lf{0) ~ lim2;_i.o Lf{x), with domain 

VL = {f: f, Lf e Ct{T}, /"(0+) = /'(0+) - /'(O-)}. 
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Figure 1: Solution of the OSP for the skew Brownian with a = l/8, /3 = l/3 and reward 
function g — {x + l)"*". Black graph corresponds to g, while the gray one corresponds to 
ipa- Va is the dotted line. 



Consider the reward function g{x) — {x + l)"*", that satisfies the RRC for xi = —1 (it 
can be seen with the same arguments considered so far) . We discuss the solution of the 
optimal stopping problem in terms of the discount factor, in particular we are interested 
in finding values of a such that the optimal threshold is the sticky point 
a different way: we fix a; = and solve the resulting equation in a. 



We use (120|) in 
We obtain 



1 



/2S y 



(O.oo) 



a{y + l)2dy 



and we find that ai — - — ^^^"^ 
ation of Theorem 13.11 



0.19 is the unique solution. It can be seen, by applic- 
that if a — ai the problem is right-sided with threshold x* — 0. 
In this case the representing measure of Va does not charge x* despite the speed meas- 
ure does. Furthermore, Theorem 4.1 can be applied to conclude that Va is differentiable 
with respect to ipa- It also should be noticed that both SF and SSF fail to hold in this 
case. This was expectable because the sufficient conditions given in [24] and [22] for the 
different types of smooth fit are not fulfilled. Another interesting thing to remark is that 
d{Va/(pa)/ds (and also d{Va/ipa)/ds) exists at 0, which is part of the conclusion of The- 
orem 2.1 in [21], despite this result is not applicable in this case. This last fact seems to 
be related with the existence of dVa/dipa 



Another approach to obtain x* = is when the strict inequality holds in ( 16 1 and also 



(191 holds. We solve (in a) equation (19) (with equality), which is 



5(0) 



w, 



:Va(o) 



i^a{y){a- L)g{y)m{dy). 



(26) 



Since the measure m{dy) has an atom at y 



from ai. Solving (26) we find the root a2 



0, the solution of the previous equation differs 
1/2. It is easy to see that for a G [ai, 012], the 



minimal x satisfying the inequality ( 16 ) is 0. Theorem 3.1 can be applied to conclude that 
the problem is right-sided with threshold x* — 0. For a G (ai,a2] we cannot apply any 
of the results of section |4j and in fact, for a S (ai,a2) none of the smooth fit principles 
is fulfilled. With a — a2 there is SF (and also SSF, since s{x) = x), but this is not 
a consequence of (]26|), it is due to the particular reward function. This example shows 



that theorems on smooth fit in section [4] only gives sufficient conditions. In table 5^ we 
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SBM: solution of the OSP depending on a 



a X* (161 (19) SF&SSF 3 dK^/rfV-a 3 d(K./v3a)/ds Fig. 



2fa: 



a € (0,Q!i) X* > = = yes yes yes 

a = ai X* = = < no yes yes 

a e {ai,a2) x* = > < no no no 

a — a2 X* = > ~ yes no no 

a e (q;2,+oo) X* < = = yes yes yes 2(e) 



summarize the information about the solution of the optimal stopping problem. We also 
give, in Figure [2| some graphics showing the solution for different values of a. 
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